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GLOSSARY 


d Spacing between dipole elements. 

Elevation ctngle of incoming or outgoing wave front, 
pgjj Azimuth angle of incoming or outgoing wave front. 

♦„ or p Relative phase along x axis. 

♦y or Py Relative phase along y axis. 

M Multiplicity or number of branches stemming from each 

node. 

Nq Nus^er of dipoles in square array. 

L Number of distribution levels when = Njj. 

°N,. N... jj Location of particular dipole. 

1 ^ 2 * L 

N Number of ^ radius in node departure angle. 

\ Order of level l<k<L. 

g Number of bits controlling phase shifts; also 

LSB prior to truncation. 

P Index integer 0, 1,2, .. . 

Np Number of phase shifts in array. 

Gjj Total wire length in single cell of multiplicity M. 

Wjj Total wire length in array of multiplicity M. 

A Expensive operational block of computer. 

B Less expensive operational block of computer. 


INTRODUCTION 


A uniqae dipole phased array antenna has been proposed* 
which utilizes micro electronics distributed throughout the 
array. The objective of this paper is to discuss the initial 
signal processing and logical organization of the distributed 
systaa. Some of the relationships between signal distribution 
methods, number of components, wire lengths, inherent errors 
and computation econ<miics are derived and their trade-offs 
discussed. Although the relationships derived are applicable 
to an ^u:ray of arbitrary size and shape, an example of a square 
4000 element eurray is used in tlie conclusion to illustrate the 
application of these relationships. 


* Harrington, J. V., R. H. Baker and J. C. James: "Digitally 

Controlled Phased Array for Reception of Siinblazer Signals and 
Solar Radar Studies," M.I.T. Center for Space Research, TR-66-3. 


1.0 POINTING EQUATIONS 

Consider the vector N pointing in the direction of wave 
propagation (i.e., normal to the plane wave front as shown in 
Figure a~ 1. When the wave front is coincident with point d 

A 

on the :c-axis, the distance fr<»a the wave front to the origin 
is N. This distance represents (in wave lengths) the relative 
phase between the point d and the origin. This phase expressed 
in radians is: 


= 2w(N) 
X 



cos 


[11 


or in terms of the azimuth angle, and elevation angle, 

p^l (see Appendix A) : 
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And the phase between a point on the y axis, d^, cmd the origin 


is: 
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The quantities 
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are referred to as phase tapers and have the advantage that they 
have no dependence on the physical quantities of the antenna 


array 
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2.0 DISTRlBOTIOl? SCHEME 

Figures lA and 2A sho:’ poscibia envisioned arrays while 
IB and 2B show these arrays rotated 45* to provide reduction 
in mathematical complexity. Envisioned signal paths* are 
shown as a series of distribut5.ori nodes. Note that the total 
electrical distance to each dipole is the same. 

If each node is considered to be in the center of a cellr 
the whole array can be thought of as a nested set of cells. It 
is convenient to define the number of subcells contained in any 
given cell as being the multiplicity, of the array. M is 
also the number or multiplicity of the branches steioming from 
each node. The total ntnaber of dipoles in the array then is 
Njj * where L is the number of distributing points or nodes 
between any dii>ole and the outside world; and, therefore, is 
also the number of different levels of cells. 

3.0 METHOD OP DISTRIBUTING PHASE CONTROL 

It has been suggested that the phase imparted to each 
dipole be controlled through the use of phase shifters located 
at each node. Thus the phase of a particular dipole would be 
the total {^ase shift accumulated by the signal as it passes 
through all nodes between the receiver and the dipole. 

L 

* E ♦, 16 ] 

^dxpole 

where k denotes the order of the node. 
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As suggested by R. H. Baker and staff 
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It becomes convenient at this point to assisme that the 
multiplicity is 4. If we assume the configuration in Figure 
IB, then each pheise shift occurs either in the x direction or 
y direction with opposing distributions from a node having 
opposite phase shifts. If we define a dipole location by the 
route a signal must take to reach it starting from the center 
of the array (first node) , then we can express directly the 
phase of that dipole as a function of its location. Let us 
further define the direction embarked fr<xa a particular node 
as N « , i.e. for N*0 along positive x axis^ N*1 along 

positive y axis, along negative x eucis, and N»3 along 
negative y axis as shown in Figure 3A. Then the location of 
a particular dipole can be given by a series of numbers each 
of which can take on values from 0 to 3. Figure 3B gives an 
example for a 64 element array. 




O-j -I n location of dipole in Figure 33 
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The relative phase at this dipole thus equals; 
Re I £ 2^“ I - 
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The above equation implies that only p and p , need to 

az el 

be specified (i.e. the pointing angles) to control the phase 
at any dipole in the array. This equation also shows that the 
phase information supplied to each node is the same except for 
a factor of 2 depending on the cell level of the node. The 
values for each node are merely the departing directions and 
are automatically hard~wired into each node. 

Thus if the two magnitude values of the expression in 
[7C] (Nj^ = even, SJj, « odd) are e::pressed as binary words and 
the phase shifters are linear, then only these two words need 
be supplied to the nodes to control the phase distribution of 
the entire array. The only calculations needed for different 
node levels are binary shifts (i.e. division by two). Even if 
the multiplicity distribution scheme (Figure lA) is not used, 
[7C] is still useful in that a distribution computor only needs 
to solve this equation to direct the array. 


4.0 QUANTIZATION 

If we were using continuous phase shifters, the relative 
phase of each dipole could be controlled exactly (within engi- 
neering tolerances). However, economical considerations indi- 
cate that quantized phase shifters might be more practical. It 
is the purpose of this section to explore some of the implica- 
tions and restrictions imposed through the use of quantized 
phase shifters. 
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Let the following binary word be a typical phase shift 
word of infinite resolution: 

0011.010111:01001 . . . 2it(3.3643 . . .) 

A A 

1 st bit qth bit 

where those bits to the left of the period denote a number of 
"2ir" phase shifts and those bits to the right denote that por- 
tion of the phase shift less tlian 2it. The denotes the 
point of guemtization. The LSB (least significant bit) above 
the guemtization shall be referred to as bit q, i.e. the LSB 
carries 2w x phase angle. 

The bits prior to the period can be ignored at each phase 
shifter in that we are working with a modular 2ir system. The 
bits following ":** will be referred to as the residue. 

If the desired pointing angle has a non- zero residue, then 
quantization will introduce a quantizing error dependent on the 
residue. The maximum peeUc error at the worst case dipole due 
to quantization is discussed in the following two sections in 
that it is akin to a tolerance in an antenna dish and also be- 
cause it can be readily calculated. 

In 4.1 axid 4.2 it is assumed that the proper phase shift 
word is supplied to each node and the quantization at one node 
does not alter the word supplied to the other nodes. 

4 . 1 TRUNCATION 

The first and perhaps easiest form of quantization is 
truncation, i.e. just ignoring the bits below the guantizaticm 
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level q. Truncation would result in the largest possible error 
at any dipole approaching 2t x x L. Noxrmalizing this error 
by the LSB, i.e. 2v x 2~^, this meucimum error approaches L 
(see^J^pendix B for detailed discussion). 

4 . 2 ROUNDING OFF 

Rounding off prior to truncation improves (reduces) the 
peak error due to quantization. The improv€snent factor approaches 
2^ for and is equal to 2 for all odd M or large even M (see 

Appendix B) . This improvement is significant enough to warremt 
rounding off prior to all truncation. 

A detailed discussion of the quantization errors due to 
rounding off is given in Appendix B. It is worthwhile noting, 
however, that the peak error for a very large array cam always 
be reduced by increasing M in that the effect of a reduction of 
nodes more than offsets the improvement factor advantage of 
M«4. 


4.3 QUANTIZING THE LOOK ANGLES 

If the pointing directions are quantized such that [8] is 
forced to hold for any integer value of P (see (7CJ , 

cos Re (e^ ^**az x 2irP (8) 

P“0, 1,2. . . 

then the residue in the phase shift words will be zero at all 
nodes and no quantization error will exist; however, a pointing 


11 - 


error will exist. If the loss in overall gain due to incorrect 
pointing is greater than it would be due to the presence of 
unfavorable residues (lack of smoothness in the phase taper 
approximation) , then a pointing quantization should not be 
employed to the full extent. 

Figure 4 shows the most favorcible (M=Njj) quantized look 

cuxgles for q*4 (i.e. phase angles quantized in — ^ (-^) incre- 

^ 16 

ments) . Note that even looking straight up# the worst pointing 

3 6 

error due to pointing quemtization is - 7 — = 2.5®. Thus even for 
a beam width as wide as 5®, this point v;ould be 3 db down. For 
M=4 there is further spread of these look angles /2 , making 
the main pointing error equal to 3 . 6 ®. 

By allowing some pointing error and some quantizing error, 
in the section 4.2 sense, an optimum distribution scheme can be 
achieved. This optimum would quantize the look angles in smaller 
increments than 2 it x 2 *”^ and allow only those residues that 
would contribute quantizing error to the last few nodes. 

Now increasing q has a very strong effect on the in^rove- 
ment of both types of errors. However, ccxnputation costs and 
distribution costs would increase with q so that q should be 
held to as small a value as is practical with respect to errors. 
The procedure for calculating the smallest allov?able value of 
q is involved but not difficult and deoends on the size of the 
array. In general, a larger array means a narrower beam width, 
and thus a larger value of q is necessary to keep the overall 
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errors within tolerance. Therefore, the computational costs 
increase at a faster rate them the array size. 


5.0 ECONOMICS 

In choosing a distribution system it becomes importcunt 
to consider the different cost aspects to arrive at some optimum 
criterion. The factors considered here are: 

1. Number of phase shifters per dipole 

2. Wire length 

3. Computation costs 


5.1 PHASE SHIFTERS PER DIPOLE 


Assuming again that the phase shifter can add both a posi 
tive phase increment and a negative phase increment, we can 
proceed to calculate the number of phase shifters per dipole 
as a factor of M and L (without this assumption multiply N 
by 2) . 

Thinking of the phase distributing method for a moment, 
it is not difficult to see that the number of phase shifters 
for M^even can be expressed as: 




N, 


- N 


“p = -f ^ ^ 2^2 

and for M=odd; 


2M 


° - 1 JL rw n 

2 M-1 
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M-1 
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M-l 


M~1 

M 
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M -1 
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A plot of for some values of M is shown in Figure 5. 
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NUMBER OF PHASE SHIFTERS (Np) 
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MULTIPLICITY 


FIG.5 


NUMBER OF PHASE SHIFTERS PER DIPOLE 
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FIG.6 

NORMALIZED WIRE LENGTH vs. MULTIPLICITY 
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Note that only for or 1^16 is there much added cost with 
regard to phase shifters using the multiplicity distribution 
appiroach . 


5.2 WIRE LENGTH 

Wire length is a real coz^sideration in that, for such a 
lazrge number of dipoles, wire length can contribute significantly 
to the overall cost. If it is assumed that all distributions 
are initiated from the geometric center of the array, then it 
is possible to express the total wire needed as a function of 
M emd Njj (see Appendix C) . The results of Appendix C are shown 
in Figure 6. Note that for = 64 a saving factor of 2^ in 
total wire length cem be realized by setting M=4 instead of 
M*64. The savings with M=4 for an arbitrarily large array is 
approximately : 


^ " »^2d(Njj-Npi) 


1 

2 


» 0.3N. 


Ill] 


For Njj = 2^2 - 4000 this saving in total wire length is nearly 
a factor of 20. 

The general expression for the wire saving as a function 
of M and Njj is given in Appendix C.3 [C-3] . 

5.3 COMPtJTER COSTS 


A cell of arbitrary M can be thought of as composed of 
a superposition of (M=4) subcells as shown in Figure 7. The 
only differences between each of the M=4 subcells are a scale 
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factor cuad a rotation. The computer block in Figure 8 shows 

how the phase word can be ccnaputed for any one of the M*4 sub- 

M M** X 

cells. Since a cell can be broken down into j ^ if M is 
odd) subcells and each subcell needs two words, the total manber 

Mjv f ) T 

of phase shift words that have to be generated is -»- ( ' — 

M In N ^ z 

for odd M) or a - '- a— words for the entire array. This means 

„ ^ " M In N 

that ^ block A's and block B's must be constructed. 

Fortunately, block A is the more expensive of the two block 
types in that it contains sin and cos generators along with scmie 
adders atnd multipliers. If the phase shifters are linear, then 
blocks 3 are of low cost in that each word is related to the 
output of block A by a small multiple (M^). If M is a power 
of 2, then the cost is effectively zero in that the outputs 
can be taken frem the same word with only a difference in the 
wiring. Non-linear phase shifters should be avoided unless the 
increased cost of block B (of which there are many) is ccaapen- 
sated for by a saving in the cost of the phase shifters. 

If a multiplicity is a power of a lower multiplicity (i.e. 
16, 64, 81, 125, etc.), then a variation in the aforementioned 
ccaaputer can result in considerable savings. Figure 9 illus- 
trates this point. Now the outputs of blocks A can be added 

to yield the phase shift words. Thus the number of A blocks 
P ”s 

(with M=Mg) becomes | — for odd) resulting in consider- 

able saving of ccanplexity and cost. The number of B blocks 

. p 

remains unchanged. If, in addition M = 4, i.e. M * 4 , then 
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only a single block A is needed with shifters replacing the other 
A blocks (Figure 10) . 

It is interesting to note that even if one were to set 
M = Njj = Mg and thus not tise the multiplicity concept in the 
distribution, the multiplicity concept is still very important 
to the efficient realization of the computer. 


6.0 SUMMARY 

There are definite advantages associated with both large 
multiplicities and small multiplicities. The advantages of 
large multiplicities are: less quantization error (smaller 

value of L), fewer phase shifters, more central control (phase 
shifts more concentrated) and the advantage of fewer connections 
and less electrical distance between a dipole and the receiver. 
On the other hand, small multiplicities have the advantage of 
less total wire length and a simpler and cheaper computer which 
in turn implies easier control. 

The final choice is dictated by the geometry and number 
of dipoles. For example, if the array is small, wire length 
ceases to be a consideration. 

6.1 LET Njj = 2^2 OR 4,000 DIPOLES 

The consideration of a particular example will illustrate 
some of the trade offs encountered in the choice of M. For 
Njj = 2 ^^, a good choice of parameters would be M=16 and L*3. 
These choices make for a simple ccxnputor (since M^ = M where 
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M *4) near min imum wire length (40 percent more than M=4 but 
s 

only about ^ the wire length for M=2^^) , only about 15 percent 
more phase shifters than optimum, and a controllable quantiza- 
tion error. I would use q=*4 (four bits to control phase shifter) 
for the last phase shift distribution (-^ maximum error, 2,000 
phase shifters); q=5 (five bits control) for the middle phase 
shifter distribution (^ meucimvim error, 128 phase shifters) and 
q>6 (six or more bit control) for the first set of 8 phase 
shifters. This distribution scheme would make the total maximum 
error due to quantization less than ^ while enabling pointing 
cuigle resolution to better than a minute of arc. 
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APPENDIX A 


DERIVATION OF RELATIONSHIPS BETWEEN p , p , AND p , p , 
(Figure A) x y az el 
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, cosp , - cosp 

“ 1 cosp 2 coip^^ (1 - sin^p^^) 

_ 1 °^^^el . cosp 

■ ^ =“P^^ (cos^p^^) + 5 eosp^^ (cos^p^j) - cosp^^ cosp^^ 


also: cosp = cosp - cos (90®-p ) = cosp , sinp 

P Sx oZ 6x d.Z 

Therefore, the relationships between p , p , p ,, and p are; 

X ^ 6 X 3.Z 

cosPj^ = cosPg^ cosp^^ [A-11 


cosPy = cosp^^ sinp^^ 


CA-2] 
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APPENDIX B 

ERPORS DUE TO QUANTIZATION 

A 9 a signal passes i^hrough a node towards a dipole, its 
phase is changed by either a positive or negative increment. 
Since signals distributed through opposite branches are changed 
by increments identical in magnitude but opposite in sign, it 
is always possible to follow a signal accumulating only positive 
increments on its way to a radiating dipole. The error at the 
worst case dipole for any arbitrary pointing angle is therefore 
just the sum of the magnitudes of the errors accumulated at each 
node. Since only that portion of the pointing angle eliminated 
by quantization (i.e. the residue) contributes to the error, 
it is sufficient to study the error as a function of this resi- 
due normalized by the quantization level. 

Figure B-IA shows the relationship between the worst case 
normalized dipole error and the normalized residue for a single 
cell (L*l) and M=4 with pure truncation. Here the error is 
merely the residue per se. Figure B-IB shows the resulting 
error when L is extended to 2. Note that this is just a super- 
position of B-lA and a double term denoting the additional 
error contributed by a second node. The remaining B-1 figures 
show what happens as L increases further. The maximum total 
error is reached when the normalized residue is just short of 
1 and in effect is equal to L as would be expected. 



FIG.B-IC L=3,M = 4 


FIG. B-ID L = 4,M = 4 
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PEAK ERROR AT V>/ORST CASE DIPOLE (due to quantization) 
WITH PURE TRUNCATION (AS A function of the normalized 
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FIG. B-2A L=l, M = 4 
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FIG. B-2 

PEAK ERROR AT WORST CASE DIPOLE (DUE to quantization) 
WITH ROUNDING OFF (prior to truncation as a function 

OF THE NORMALIZED RESIDUE) 
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Pigxire B2-A shows the effect on the error as a result of 
rounding off prior to truncation for L=l. As expected, the 
peak value is Figure B2-B shows the maximum error along 
either the x-axis or y-axis for L*2. As in the case of Figure 
Bi-B, this figure is a superposition of Figure B2-A emd a double 
error term resulting from the second node. This double term 
does not peak where the first error does. The peak error is 
not 2 X ^ 1, therefore, but rather ^ “ "J* Figures B2-C, 

etc. show this error for larger values of L. It is interesting 
to note that, for larger values of L, the total error approaches 
X L and not x L as might be expected. The improvement then 
along the axes resulting from rovinding off prior to truncation 
is a factor of 3. 

The off-axis errors are larger than the on-axis errors 
in that the x-cixis residues and y-axis residues are uncorrelated. 
The worst possible case is generated by alternating axes, i.e. 
as close to the 45® diagonal as the discreteness of the array 
will allow. For large L and M=4, the diagonal error approaches 
0.4 X L which is a factor of improvement over pure trunca- 
tion. 

For M=9 the worst case peak errors occur along the eucis. 
The total error is again found by a superposition of the errors 
contributed at each node. Since all the contributions peak for 
the same residue, this is the maximum peak error which an off- 
axis cannot exceed. For M=25, 49 . . . the maximum peak errors 
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will also occur along the axis and each node will contribute 
eui error of 0.5 for a total meucimum peak error of 0.5L. 

For M=16 very little cross corirelation remains between 
the error terms , and thus the diagonal meucimum peak error is 
almost 0.5L (actual 0.47L). For M=36, 64, 100 . . . the mcuci- 
mum peak errors are also along the diagonals and are effectively 
0.5L. 
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APPENDIX C 


CALCULATION OF WIRE LENGTHS VERSUS M AND N^ 


Let us define the wire length in the small cell of multi** 
plicity M as Gj^d (where d is the dipole spacing and G^j is a 
geometric constant depending on M) . Then the contribution to 
the total wire length due to all of the small cells is — — 
and the total wire length can be expressed in a series where 
each term is the contribution due to those cells at a particular 
level. 


W, 
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1 + 


1 1 
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+ 
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w, 
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M-M^ 


lC-2] 


or normalizing by all terms which do not depend on the distri- 
bution method: 
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C,1 FOR M = Np FOR LARGE ARRAYS 


An approximation to the total wire length with M«Njj for 
large arrays can be calculated by multiplying the densities of 
dipoles , p, in a delta area weighted by the wire length to the 
middle of the array and then integrating over the entire array. 
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[C-4] 


or 


W, 


N, 


N, 


» 0.382 dN, 


[C-5] 


The significance of this result is that W„ increases as 
where, if M is held to some constant value, increases 
linearly with Njj. Thus for large arrays a considerable savings 
in V7ire length and therefore cost can be accomplished by holding 
M to scxne constant value (i.e. 4, 9, etc.) 
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C.2 W„ FOR ARBITRT^Y VALUES OP M 

It is of interest to compare [C-5] with tC-2] to see how 

large the array must be to make [C-5] valid. Exact calculated 

1 

values for for M=4, 9 . . . are plotted against M in 

Figure 6, page 14, as a series of dots. as calculated 

by combining [C-5] and [C-3] is shown in Figure 6 as a series 
of X's and are given in lC-6] . 

— 1 = 0.382 “T [C-6] 

M^-l 


The approximation is quite good for all but the lowest 
values of H. It is worthwhile then to substitute [C-6] into 
[C-3] to obtain an equation for V1^ dependent entirely on known 
quantities. The resulting approximation is within two percent 
for M greater than 9 and within ten percent for M equal to 4 
or 9. 



(0.382) 


dM(Njj-Njj^ 

(M^-1) 


C.3 SAVING IN WIRE LENGTH FOR M ji^ Nj^ 


[C-7] 


One final useful equation is a closed form expression 
for the saving in wire length as a function of M compared to 
M = Njj. This equation is found by dividing [C-5] by [C-7] : 




<V-» _ "d (mU) 

^ M “ (Sq^-I) 

(m2-1) 
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for large M 
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